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ABSTRACT 

We use the 4 Ms Chandra Deep Field-South (CDF-S) survey to place direct constraints on the ubiquity of 
z«2 heavily obscured AGNs in K < 22 BzK-selected galaxies. Forty seven (w 21%) of the 222 Bz/T-selected 
galaxies in the central region of the CDF-S are detected at X-ray energies, 1 1 (w 5%) of which have hard 
X-ray spectral slopes (r £ 1) indicating the presence of heavily obscured AGN activity (A/h ^ 3 x 10 23 cm -2 ). 
The other 36 X-ray detected BzK galaxies appear to be relatively unobscured AGNs and starburst galaxies; we 
use X-ray variability analyses over a rest-frame baseline of w 3 years to further confirm the presence of AGN 
activity in many of these systems. The majority (seven out of 1 1) of the heavily obscured AGNs have excess 
infrared emission over that expected from star formation (termed "infrared-excess galaxies"). However, we find 
that X-ray detected heavily obscured AGNs only comprise » 25% of the infrared-excess galaxy population, 
which is otherwise composed of relatively unobscured AGNs and starburst galaxies. We find that the typical 
X-ray spectrum of the heavily obscured AGNs is better characterized by a pure reflection model than an ab- 
sorbed power-law model, suggesting extreme Compton-thick absorption (A/h ^ 10 24 cm -2 ) in some systems. 
We verify this result by producing a composite rest-frame 2-20 keV spectrum, which has a similar shape as a 
reflection-dominated X-ray spectrum and reveals an emission feature at rest-frame energy » 6.4 keV, likely to 
be due to Fe K. These heavily obscured AGNs are likely to be the distant analogs of the reflection-dominated 
AGNs recently identified at z ~ with > lOkeV observatories. On the basis of these analyses we estimate 
the space density for typical (intrinsic X-ray luminosities of L2-i0keV ~ 10 43 erg s -1 ) heavily obscured and 
Compton-thick AGNs at z ~ 2. Our space-density constraints are conservative lower limits but they are already 
consistent with the range of predictions from X-ray background models. 

Subject headings: galaxies: active — galaxies: high-redshift — infrared: galaxies — X-rays: galaxies — 
ultraviolet: galaxies 



1. INTRODUCTION 

Deep X-ray surveys have provided a penetrating probe of 
Active Galactic Nuclei (AGN) out to z ~ 5 (e.g., Brandt & 
Hasinger 2005; Silverman et al. 2008; Brusa et al. 2009; 
Brandt & Alexander 2010), identifying obscured and unob- 
scured AGN activity in a modest fraction of the field-galaxy 
population ( > 5-10%; e.g., Lehmer etal.2005, 2008; Xue 
etal. 2010). However, there is overwhelming evidence that 
a large fraction of the heavily obscured AGN population 
(A/ H £ 3 x 10 23 cm- 2 ) remains undetected in even the deep- 
est X-ray surveys (e.g., Worsley et al. 2005; Hickox & Marke- 
vitch 2006; Tozzi etal. 2006; Treister etal. 2006; see §1 of 
Alexander et al. 2008). Distant heavily obscured AGNs are 
predicted by many theoretical models and simulations to rep- 
resent an important phase in the evolution of distant dust- 



enshrouded galaxies, where the rapidly growing central su- 
permassive black hole (BH) is hidden from view (e.g., Fabian 
1999; Granato etal. 2006; Hopkins et al. 2006). Therefore the 
identification of the most heavily obscured AGNs could be 
more than just a book-keeping exercise - without having ob- 
servations sensitive to their identification we may miss a cru- 
cial BH growth phase. 

Weak (faint or undetected) X-ray emission from luminous 
AGNs is likely to be due to the presence of large amounts of 
dust/gas, sometimes exceeding A/h ^ 10 24 cm (i.e., Comp- 
ton thick; e.g., Matt etal. 2000; Comastri 2004; Delia Ceca 
et al. 2008; Murphy & Yaqoob 2009). Strong support for 
this statement comes from the tight correlation between the 
optical and X-ray emission of unobscured quasars (e.g., Vig- 
nali, Brandt, & Schneider 2003; Steffen et al. 2006; Gibson 
et al. 2008), which suggests that all luminous AGNs are intrin- 
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sically bright at X-ray energies. Although weak at X-ray ener- 
gies, these heavily obscured AGNs should still be detected in 
deep mid-to-far infrared (IR; rest-frame wavelength > 2 /urn) 
observations due to the presence of dust heated by the hid- 
den AGN. Indeed, a number of studies have revealed large 
populations of X-ray undetected IR-bright galaxies at z « 2, 
which may host heavily obscured, potentially Compton thick, 
AGN activity (e.g., Daddi et al. 2007a; Fiore et al. 2008, 2009; 
Georgantopoulos etal.2008, 2011; Treister etal.2009; Geor- 
gakakis et al. 2010). Many of these studies have employed 
X-ray stacking analyses of X-ray undetected IR galaxies to 
identify hidden AGN populations statistically, where the de- 
tection of a hard X-ray spectral slope (Y ^ 1) in the stacked 
data provides compelling evidence for the presence of heav- 
ily obscured AGN activity in at least a fraction of the stacked 
sources. 1 Under the assumption that all of these X-ray unde- 
tected IR galaxies host Compton-thick AGN activity, the im- 
plied space density of these hidden AGNs would exceed those 
of the Compton-thin AGN population by a factor of £ 2, im- 
plying that a much larger fraction of distant luminous AGNs 
are Compton thick than found locally ( <, 25-50% of local 
AGNs appear to be Compton thick; see Table 2 in Burlon et al. 
2010; Risaliti et al. 1999; Guainazzi et al. 2005). These dis- 
coveries provide some support for the hypothesis that the ma- 
jority of distant BH growth was more heavily obscured than 
that found locally (e.g., La Franca et al. 2005; Treister & Urry 
2006; Hasinger 2008). 

However, results from X-ray stacking analyses need to be 
treated with caution since they only provide an average sig- 
nal, leaving significant uncertainties about the overall distri- 
bution of source properties. Therefore, before strong conclu- 
sions can be derived from these studies, at least two key ques- 
tions need to be addressed: (1) what fraction of the X-ray 
stacked signal is "contaminated" by star-forming galaxies, 
which are found to comprise at least a fraction of the can- 
didate heavily obscured AGN population (e.g., Donley et al. 
2008; Murphy et al. 2009; Fadda et al. 2010)? (2) what frac- 
tion of the heavily obscured AGNs are absorbed by Comp- 
ton thick (as opposed to Compton thin) material? The most 
direct way to address these questions is with deeper X-ray 
data, which will (1) reveal the X-ray properties of individ- 
ual IR-bright galaxies that were previously contributing to the 
stacked X-ray signal, (2) allow for more detailed X-ray spec- 
tral investigations of the X-ray detected IR-bright galaxies to 
search for the signatures of heavily obscured and Compton- 
thick absorption (e.g., the identification of a strong reflection 
component at rest-frame > 10 keV; the detection of a high 
equivalent width Fe K emission line; e.g., Matt etal. 1996, 
2000; Tozzi et al. 2006; Georgantopoulos et al. 2009; Co- 
mastri et al. 201 1 ; Feruglio et al. 201 1 ; see Murphy & Yaqoob 
2009), 2 and (3) improve stacking constraints of X-ray un- 
detected populations. The presence of large amounts of ab- 
sorption can also be indirectly inferred from the identifica- 
tion of luminous AGN emission lines and an IR-emitting hot- 
dust AGN continuum in X-ray weak systems (e.g., Bassani 

1 The spectral slopes of X-ray emission from star-formation processes are 
typically T > 1 (e.g., Kim etal. 1992; Ptak etal. 1999; Berghea etal. 2008; 
Iwasawa etal. 2009). 

2 We note that high signal-to-noise ratio data is required to accurately dis- 
tinguish between the X-ray spectral properties of Compton-thin AGNs with 
Nr ~ (5-10) x 1() 23 cm~ 2 and Compton-thick AGNs with Nr > 10 24 cm" 2 
(e.g., Murphy & Yaqoob 2009; Yaqoob etal. 2010). Therefore our adopted 
definition of an Compton-thick AGN is not comprehensive but it is conven- 
tional; see the MYTorus manual at www.mytorus.com for a detailed review. 



etal. 1999; Krabbe etal. 2001; Lutz etal. 2004; Alexander 
etal. 2005b, 2008; Heckman etal. 2005; Gandhi etal.2009; 
Bauer et al. 20 1 0; Donley et al. 20 1 0; Gilli et al. 20 1 0; Vignali 
etal. 2010; Goulding etal. 2011). Greater reliability in the 
identification of heavily obscured and Compton-thick AGNs 
is made when considering multiple diagnostics that cross- 
check each other, particularly those that probe different AGN 
regions. 

In this paper we use the deepest X-ray observations avail- 
able (the 4 Ms Chandra exposure of the Chandra Deep 
Field-South; CDF-S; Xue etal. 2011) to extend the analy- 
ses of Daddi et al. (2007a), which employed X-ray stack- 
ing techniques to study X-ray undetected IR-bright z ~ 2 
galaxies in the shallower 1 Ms CDF-S observations. Daddi 
etal. (2007b) utilized the BzK photometric-selection tech- 
nique (Daddi etal. 2004) to identify K < 22 galaxies at 
z ~ 1.4-2.6 and classified objects based on the ratio of 
mid-IR (24 jum) to extinction-corrected ultra-violet (UV; 
rest-frame 1500 A) star-formation rates (SFRs). The BzK 
photometric-selection technique provides an effective iden- 
tification of massive galaxies (w 10 10 — 10 n M ; see Daddi 
etal. 2007b; McCracken et al. 2010). BzK galaxies with a sig- 
nificant excess of IR emission over that predicted from the 
extinction-corrected UV SFRs were classified as "IR-excess' 
galaxies" [log(SFR(mid-IR+UV)/SFR(UV,corr))> 0.5] while 
BzK galaxies with comparable mid-IR and extinction- 
corrected UV SFRs were classified as "IR-normal galax- 
ies" [log(SFR(mid-IR+UV)/SFR(UV,corr))< 0.5]. From 
stacking the X-ray data of the X-ray undetected galaxies, 
Daddi et al. (2007a) obtained distinctly different X-ray spec- 
tral slopes for the IR-excess galaxies (r w 0.9) and IR- 
normal galaxies (r m 1.8). The flat X-ray spectral slope 
found for the IR-excess galaxies indicates that a fraction of 
the X-ray undetected IR-excess galaxy population host heav- 
ily obscured AGN activity, some of which may be Comp- 
ton thick; by comparison, the stacked X-ray emission from 
the X-ray undetected IR-normal galaxies is consistent with 
that expected from star formation. Under the assumption 
that all of the IR-excess galaxies are Compton-thick AGNs, 
Daddi et al. (2007a) estimated a space density of <t> « 2.6 x 
10" 4 Mpc~ 3 for z « 2 Compton-thick AGNs with L 2 -10keV > 
1 42 — 1 43 erg s _1 . With the deeper Chandra data from 
Xue etal. (2011) we can now better characterize the X-ray 
properties of both the X-ray detected and X-ray undetected 
BzK galaxies and improve constraints on the ubiquity of dis- 
tant heavily obscured and Compton-thick AGNs. We adopt 
H Q = 71 km s~ l Mpc-\ Q. M = 0.27, and il A = 0.73 through- 
out. The Galactic absorption toward the CDF-S region is 
N H = 8.8 xlO 19 cm" 2 (Stark et al. 1992). All given mag- 
nitudes are based on the Vega-magnitude system. 

2. DATA AND STACKING PROCEDURES 
2.1. Galaxy sample 

We use an updated version of the z ~ 2 galaxy samples 
generated by Daddi etal. (2007b) in the CDF-S field. Due 
to small refinements in the optical-mid-IR photometry and 
revised redshift estimates, our updated sample is slightly dif- 
ferent to that used in Daddi etal. (2007 a,b). However, qual- 
itatively, these samples are the same as those used in Daddi 
etal. (2007 a,b) and have the same global properties. 

To provide a good compromise between excellent X-ray 
sensitivity and a large number of galaxies, we have only ex- 
plored the X-ray properties of BzK galaxies that lie within 5. '5 
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of the average Chandra aimpoint. Within this region there are 

76 objects classified as IR-excess galaxies and 146 objects 1000 

classified as IR-normal galaxies. 



2.2. X-ray matching 

We matched the parent BzK galaxy sample of 222 objects 
to the 4 Ms Chandra catalogs of Xue et al. (201 1) using a 1"5 
search radius. In total, 47 BzK galaxies have an X-ray coun- 
terpart, 28 of which are classified as IR-excess galaxies and 
19 of which are classified as IR-normal galaxies. The median 
X-ray-,K-band position offset is 0"4, which agrees with the 
median uncertainty of the X-ray source positions (which cor- 
respond to the 68% confidence level; see Fig. 6a and Eqn. 2 in 
Xue et al. 201 1); on the basis of our matching parameters we 
expect w 1.6 spurious matches. Twenty of these BzK galax- 
ies were detected in the 1 Ms CDF-S catalogs of Alexander 
et al. (2003) and a further 27 are now detected in the 4 Ms 
CDF-S catalogs of Xue etal. (2011). The properties of the 
X-ray detected BzK galaxies are presented in Table 1 . 

Thirty of the X-ray detected BzK galaxies have spectro- 
scopic redshifts: 27 from optical spectroscopy (predomi- 
nantly from VLT observations with the FORS1, FORS2, and 
VIMOS instruments; Appenzeller et al. 1998; Le Fevre 
et al. 2003) and 3 from Spitzer-IRS mid-IR spectroscopy; 
see Table 1 for the spectroscopic-redshift references. The 
median spectroscopic redshift is z S p ec = 1.78±0.35. 3 The 
other 17 X-ray detected BzK galaxies have photometric red- 
shifts and a median redshift of z p hoto = 2.18 ± 0.50; see Ta- 
ble 1 for the photometric-redshift references. The median 
absolute uncertainty between the photometric and spectro- 
scopic redshifts for the X-ray sources with spectroscopic red- 
shifts is |Az|/(l +z sp ec) = 0.02, where Az = z photo - z spec ; we 
get comparable results if we use the Cardamone et al. (2010) 
photometric-redshift catalog (|Az|/ (1 +z S pec) = 0.02). 

2.3. X-ray spectroscopy 

We extracted and analysed the X-ray spectra of the X-ray 
detected BzK galaxies to provide greater insight into their 
intrinsic properties. The X-ray spectra were extracted us- 
ing ACIS_EXTRACT (Broos et al. 2010) as part of the X-ray 
catalog construction in Xue etal. (2011); see §3.2 of Xue 
et al. (201 1) for more details. 

Due to the limited counting statistics for the majority of the 
heavily obscured AGNs ( 80 counts in the 0.5-8 keV band; 
see Table 2), which are the primary focus of this paper, the 
X-ray spectral analyses were predominantly performed using 
the C-statistic (Cash 1979). The C-statistic is calculated on 
the unbinned data and is therefore ideally suited to low-count 
sources (e.g., Nousek & Shue 1989). However, to provide 
consistency checks on these results, we also performed X-ray 
spectral analyses of the brightest X-ray sources (> 200 counts 
in the 0.5-8 keV band) using % 2 statistics; in these analyses 
we grouped the X-ray data into 20 counts per bin. All fit pa- 
rameter uncertainties are quoted at the 90% confidence level 
(Avni 1976). 

2.4. X-ray variability 

We analysed the X-ray variability of the X-ray detected BzK 
galaxies to look for nuclear activity. The CDF-S observa- 
tions were split into four epochs, each approximately 1 Ms 

3 The error on the median is the median absolute deviation (MAD), which 
is a robust estimator of the spread of the sample: MAD = 1.48 x median(|jr — 
median(jt) | ) ; see § 1 .2 of Maronna et al. (2006). 
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FIG. 1. — Ratio of star-formation rates (SFRs; mid-IR and extinction- 
corrected UV) for the z ~ 2 galaxies studied by Daddi et al. (2007a,b) ver- 
sus rest-frame 8 fim luminosity. Triangles correspond to X-ray detected BzK 
galaxies, crosses indicate the BzK galaxies with hard X-ray spectral slopes 
(heavily obscured AGNs; see Fig. 2), open squares correspond to the X-ray 
undetected BzK galaxies, and the open circles correspond to the spectroscop- 
ically identified (optical and mid-IR wavelengths) Compton-thick AGNs at 
z ~ 2 from Alexander et al. (2008) in the 2 Ms Chandra Deep Field-North 
survey (Alexander etal. 2003). The threshold between IR excess and IR nor- 
mal galaxies defined by Daddi et al. (2007a,b) is indicated by the dashed line. 



long: 2000, 2007, 2010a (March - May), and 2010b (May - 
July). Within each epoch, the observations were merged and 
photometry was measured using ACIS_EXTRACT (Broos et al. 
2010), as described in detail in Xue et al. (2011). 

We apply a % 2 test to determine if a source is variable by 
comparing the variability observed between observations to 
that expected from Poisson statistics. The test statistic will 
follow a % 2 distribution except at low count rates, where the 
errors are larger than expected from a Gaussian distribution. 
In the low-count regime, the test statistic is smaller than ex- 
pected and does not follow the % 2 distribution. We construct 
a Monte Carlo simulation to determine what distribution the 
test statistic should follow for each source, following the pro- 
cedure of Paolillo et al. (2004). The observed test statistic is 
compared to the simulated distribution to determine the prob- 
ability (Pjo) that the observed variability is due to Poisson 

noise. The % 2 and P^i values are listed in Table 2. A source 
is considered variable if it has > 20 counts in the 0.5-8 keV 
band and P % 2 < 5%. 

The normalized excess variance (a 2 nm ; Nandra et al. 1997) 
measures how strongly each source varies in excess of mea- 
surement error. Since it is more likely that a variable galaxy 
harbors an AGN if the observed variability exceeds that ex- 
pected from a population of X-ray binaries, we estimate the 
amount of variability expected from high-mass X-ray bina- 
ries (HMXBs). HMXB variability can be estimated from 
a galaxy's star formation rate (see Fig. 8 of Gilfanov et al. 
2004). To use this relation we adopt a 2 ms = 0.09 as the typical 
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variability for an individual X-ray binary (equivalent to 30% 
fractional rms; see §6 of Gilfanov 2010). This value is the 
maximum variability expected for a HMXB population. Vari- 
ability strength will increase with the length of the timescale 
over which it is measured (e.g., Nandra et al. 1997), so adopt- 
ing the maximum HMXB variability is appropriate for CDF-S 
sources at z ~ 2, where rest-frame timescales are w 3 years. 
The SFR of the BzK galaxies was estimated from their UV 
and IR luminosities (Xue et al. 2010). Applying the 0^ im - 
SFR relations (Gilfanov et al. 2004), we calculate an upper 
limit to the HMXB contribution to the variability, reported as 
^hmxb m Table 2. A source is considered to be an X-ray vari- 
able AGN if it is found to be variable, as prescribed above, and 
also has excess variability over that expected from HMXBs. 
On the basis of these criteria, 13 of the 27 sources with > 20 
counts in the 0.5-8 keV band are found to be variable AGN. 

2.5. X-ray stacking analyses 

We used X-ray stacking analyses to constrain the aver- 
age X-ray properties of the X-ray undetected BzK galaxy 
populations. In our X-ray stacking analyses we adopted 
the procedure of Lehmer et al. (2008), which takes a differ- 
ent approach from the Worsley etal. (2005) method used by 
Daddi et al. (2007a). Both procedures stack the X-ray data 
of the selected sources but Lehmer et al. (2008) determine the 
background counts using large source-free apertures local to 
each source while Worsley et al. (2005) determine the back- 
ground counts from a large number of randomly placed aper- 
tures around the source (i.e., a Monte-Carlo approach). We 
tested both procedures on our datasets and achieved statisti- 
cally consistent results. The major advantage of the Lehmer 
etal. (2008) approach over that of Worsley etal. (2005) is 
computational speed. In the stacking analyses, we used a fixed 
aperture of 1 "5 radius and determined background counts in 
25" x 25" source-free regions local to each source. We ap- 
plied aperture corrections to the net stacked count rates fol- 
lowing §4.2 in Lehmer et al. (2008). 

A general concern in stacking analyses is that a few sources 
can dominate the stacked signal. To guard against this, we 
randomly selected 80% of the objects in each stacking analy- 
sis sample and stacked their properties. For each sample we 
performed this procedure 10,000 times, randomly selecting 
80% of the objects for each iteration, to generate a distri- 
bution of the stacked properties. We found that the overall 
properties obtained from the stacking analyses were in good 
agreement with the overall properties found from the stacking 
analysis trials, indicating that bright sources do not dominate 
the stacked signal. 

3. ANALYSES AND RESULTS 

We focus our analyses toward (1) characterizing the X-ray 
spectral and variability properties of X-ray detected BzK 
galaxies to identify the presence of AGN activity (see §3.1), 
and (2) performing X-ray stacking analyses of the X-ray un- 
detected BzK galaxies in the 4 Ms CDF-S observations (see 
§3.2). With the results of these investigations we re-evaluate 
estimates of the space density of distant heavily obscured and 
Compton-thick AGNs (see §3.3). 

3.1. X-ray detected BzK galaxies 

In Fig. 1 we compare the properties of the X-ray detected 
BzK galaxies to the overall BzK galaxy population. The 
X-ray detected BzK galaxies cover a wide range in SFR ratio 
(log(SFR(mid-IR+UV)/SFR(UV,corr)) and rest-frame 8 fim 
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FIG. 2. — Rest-frame 2-10 keV luminosity versus X-ray spectral slope (r is 
in the observed-frame 0.5-8 keV, typically rest-frame 1 .5-24 keV). The prop- 
erties of z, ~ 2 SMGs hosting AGN activity (open pentagons; from Alexan- 
der et al. 2005a) and X-ray detected Bz,K galaxies (filled triangles; identi- 
fied here) are compared to well-studied local starbursts (open stars) and lo- 
cal Compton-thick AGNs (open squares). The properties of the local star- 
burst galaxies and Compton-thick AGNs are taken from Ranalli et al. (2003), 
Matt etal. (1997), and Vignati et al. (1999). The vertical dashed lines indicate 
the difference between the observed X-ray luminosity and the absorption- 
corrected luminosity for the SMGs and Compton-thick AGNs (from Matt 
etal. 1997; Vignati etal. 1999; Alexander etal. 2005a; Feruglio et al. 2011); 
the intrinsic X-ray luminosity constraint for J033235. 7-274916 is from Fer- 
uglio et al. (201 1), which is one of the heavily obscured AGNs in our sample. 
The range of X-ray spectral slopes found for typical AGNs and star-forming 
galaxies are illustrated (horizontal solid lines; constraints derived from Kim 
etal. 1992; Nandra & Pounds 1994; Maiolino etal. 1998; Ptak etal. 1999; 
Berghea et al. 2008). Eleven of the X-ray detected BzK galaxies have flat 
X-ray spectral slopes (F^, 1) and are classifed as heavily obscured AGNs 
(crosses). 



luminosity. However, when compared to the X-ray unde- 
tected BzK galaxies, the X-ray detected systems have char- 
acteristically higher median SFR ratios and rest-frame 8 fim 
luminosities (X-ray detected: SFR ratio of m 4.5 ± 4.9 and 
log(L 8 ^m/L Q ) w 11.1 ±0.4; X-ray undetected: SFR ratio 

of « 1.4 ± 1.4 and log(L 8 pm/L @ ) ~ 10.6 ± 0.4). 4 Over- 
all, w 37% of the IR-excess galaxy population and ps 13% 
of the IR-normal galaxy population are now detected in the 
4 Ms Chandra exposure, indicating a close connection be- 
tween the production of X-ray emission and the presence of 
excess (or luminous) IR emission. Indeed, the X-ray detected 
fraction rises as a function of rest-frame 8 fum luminosity: 
w 3%, w 14%, and « 51% of the log(L 8 m /L & ) = 10.0-10.5, 
log(L 8(Um /L Q ) = 10.5-11.0, and log(L 8 A , m /L ) > 11.0 sys- 
tems are detected at X-ray energies, respectively. Given these 
results, X-ray observations an order of magnitude deeper than 
those obtained here (only likely to be attainable with the next 
generation of X-ray observatories; e.g., Generation-X; Wolk 
et al. 2008) are required to individually detect X-ray emission 



4 Lgu 



10 L is comparable to Lir J 10 i0 for the spectral energy 



distribution corrections adopted in Daddi et al. (2007b). 
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FIG. 3. — Rest-frame 2-10 keV luminosity versus rest-frame 8 /jm luminos- 
ity. The symbols have the same meaning as in Figs. 1 & 2; in addition, the 
small filled circles show the X-ray detected BzK galaxies with Spitzer-IRS 
spectroscopy. The X-ray-8 jum luminosity ratio for local starburst galaxies 
(long-dashed line) is taken from the X-ray-12 /jm luminosity ratio of Krabbe 
etal. (2001) and converted to 8 fim assuming the M 82 spectral energy distri- 
bution. The intrinsic X-ray-8 /jm luminosity ratio for local AGNs (solid line) 
is taken from Lutz etal. (2004) and converted to 8 fjm, assuming the AGN- 
dominated galaxy NGC 1068; the dotted line indicates the observed X-ray- 
8 fira luminosity ratio predicted for Compton-thick AGNs. The rest-frame 
8 /jm luminosity for the X-ray detected BzK galaxies is calculated from the 
24 /jm flux density, with small ^-corrections applied (see Daddi etal. 2007a), 
while the rest-frame 8 fim luminosity for the local starburst galaxies is cal- 
culated using the mid-IR spectroscopy of Rigopoulou etal. (1999) and Lutz 
etal. (2003). The vertical dashed lines indicate the difference between the 
observed and intrinsic X-ray luminosity for the z ~ 2 Compton-thick AGNs 
from Alexander et al. (2008) and Feruglio et al. (2011); the intrinsic X-ray 
luminosity constraint for J033235. 7-2749 16 is from Feruglio et al. (2011), 
which is one of the heavily obscured AGNs in our sample. 



from the majority of the lowest-luminosity systems. 

3.1.1. Classification of the X-ray emission 

In Fig. 2 we plot the rest-frame 2-10 keV luminosity ver- 
sus X-ray spectral slope of the X-ray detected BzK galax- 
ies and compare them to well-studied local starburst galax- 
ies, Compton-thick AGNs, and z«2 submillimeter-emitting 
galaxies (SMGs) hosting AGN activity; the rest-frame 
2-10 keV luminosities are calculated from the observed- 
frame 0.5-2 keV fluxes assuming r = 1.8 for the small K 
corrections. Eleven (« 23%) of the 47 X-ray detected BzK 
galaxies have flat X-ray spectral slopes with T 1 (eight 
have r i$j 0.5) and are classified here as "heavily obscured 
AGNs". On the basis of the X-ray properties of the z ~ 2 
SMGs hosting AGN activity (Alexander et al. 2005a), the flat 
X-ray spectral slopes for these heavily obscured AGNs sug- 
gest absorbing column densities of Nh ~3 x 10 23 cm~ 2 and 
some may be Compton thick; see Fig. 2. The absorption cor- 
rections for such heavily obscured AGNs in the rest-frame 
2-10 keV band are large and would imply intrinsic luminosi- 
ties of Z^-iOkeV ~ 10 erg s -1 . Six of these heavily obscured 
AGNs are not detected in the 1 Ms CDF-S data used in Daddi 



etal. (2007a). 

On the basis of the X-ray luminosity and X-ray spectral 
slope, the X-ray emission from the majority of the other 36 
X-ray detected BzK galaxies is likely to be due to either rel- 



atively unobscured AGN activity (Nh 10 —10 cm ) or 
star formation. Nine of these objects have L\ > 10 43 erg s -1 , 
five of which are identified as AGNs based on optical spec- 
troscopy (see Table 1), and are classified here as "luminous 
AGNs". The other 27 X-ray detected BzK galaxies are likely 
to be lower-luminosity AGNs or X-ray luminous starbursts 
and are classified here as "low-luminosity X-ray systems"; as 
shown in Fig. 3 and §3.1.2, the X-ray-IR luminosity ratios of 
these systems are also similar to those expected for starburst 
galaxies or z ~ 2 AGNs with L x ~ 10 42 — 10 43 erg s _1 (e.g., 
Krabbe etal. 2001; Mullaney etal. 2010; Shao etal. 2010). 
Twenty one of these low-luminosity X-ray systems are not 
detected in the 1 Ms CDF-S data used in Daddi etal. (2007a). 

Four of the 27 low-luminosity X-ray systems are detected 
in both the 2-8 keV and 0.5-2 keV bands and their X-ray 
spectral slopes are F w 1.2-2.0. The other 23 low-luminosity 
X-ray systems are undetected at 2-8 keV and we can only 
provide accurate constraints on their X-ray spectral slopes us- 
ing X-ray stacking analyses. Stacking the X-ray data for these 
23 systems, following the procedure outlined in §2.5, we ob- 
tain significant detections in both the 2-8 keV (S/N = 8.3) 
and 0.5-2 keV bands (S/N = 29.7), which corresponds to an 
average X-ray spectral slope of F = 1.5 ± 0.1. The compar- 
atively steep X-ray spectral slope from this population could 
be due to either star formation or AGN activity and does not 
provide significant new insight into the composition of the 
low-luminosity X-ray systems. 

We can further characterize the X-ray detected BzK galax- 
ies using X-ray variability analyses. The identification of 
significant X-ray variability over that expected from star- 
formation processes will indicate the presence of an AGN; 
see §2.4. Overall, we find that 13 (« 48%) of the 27 X-ray de- 
tected BzK galaxies with reasonable-quality X-ray data (> 20 
X-ray counts) show excess variability over that expected from 
star-formation processes; the sources show variability by fac- 
tors of 1.4-4.3 (see Table 2). Nine of these 13 systems 
had already been classified as AGNs: six are luminous AGNs 
and four are heavily obscured AGNs. However, three of the 
variable sources are classified as low-luminosity X-ray sys- 
tems (« 40% of those with reasonable-quality X-ray data), 
unambiguously identifying the presence of AGNs in at least a 
fraction of the low-luminosity X-ray systems. 

3.1.2. X-ray-infrared properties 

In Fig. 3 we plot the rest-frame 2-10 keV luminosity versus 
rest-frame 8 pm luminosity of the X-ray detected BzK galax- 
ies and compare them to well-studied local starburst galaxies 
and Compton-thick AGNs. This figure can help character- 
ize the X-ray-8 pm luminosity ratio and provide constraints 
on the intrinsic luminosity of the heavily obscured AGNs. For 
example, under the assumption that the 8 pm emission is dom- 
inated by AGN activity, the X-ray-8 pm luminosity ratios of 
the heavily obscured AGNs suggest that they have intrinsic 
X-ray luminosities of L2-I0kev ~ 3 x 10 43 — 10 45 erg s (see 
Fig. 3). The presence of absorption at 8 pm, as predicted 
by radiative-transfer modeling of clumpy AGN obscuration 
(e.g., Nenkova et al. 2008), will increase these instrinsic X- 
ray luminosity estimates; however, see Lutz et al. (2004) and 
Gandhi et al. (2009) for observational constraints suggesting 
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that obscured AGNs do not typically suffer significant nuclear 
absorption at infared wavelengths. Under the assumption that 
the 8 /jm emission is AGN dominated, the X-ray-8 /am lumi- 
nosity ratios for all of the heavily obscured AGNs, except for 
J033222. 5-274603, are also consistent with those expected 
for Compton-thick AGNs. 

However, the interpretation of the X-ray-8 /jm luminosity 
ratio is complicated by the absence of mid-IR spectroscopy 
for the majority of the X-ray detected BzK galaxies, which 
would directly measure the contributions from star formation 
and AGN activity at rest-frame 8 /urn. Nine of the X-ray 
detected BzK galaxies have Spitzer-IRS spectroscopy, three 
of which are classified as X-ray luminous AGNs and six of 
which are classified as low-luminosity X-ray sources; see Ta- 
ble 1 and Fig. 3. With the exception of the luminous AGN 
J033237.7-275212 (which is AGN-dominated at 8 /urn; Don- 
ley et al. 2010), all of the other 8 systems are star-formation 
dominated at 8 /urn, including two of the luminous AGNs 
(Teplitz et al. 2007; Fadda et al. 2010). If the heav- 
ily obscured AGNs are the absorbed counterparts of the lu- 
minous AGNs, as suggested by the X-ray spectral analyses 
(see §3.1.3), then we would expect many of them to also be 
star-formation dominated at 8 /xm and, therefore, the range 
of intrinsic X-ray luminosities estimated above (Z-2-iOkeV « 
3 x 10 43 -10 45 erg s _1 ) are upper limits; see §3.1.4 for further 
estimates of the intrinsic X-ray luminosities. 

The majority of the heavily obscured AGNs are identified as 
IR-excess galaxies (seven of the 1 1 systems); see Fig. 1 . How- 
ever, the majority of the X-ray detected IR-excess galaxies are 
not heavily obscured AGNs: only 7 (« 25%) of the 28 X-ray 
detected systems are heavily obscured AGNs, 7 (« 25%) are 
luminous AGNs, and 14 (« 50%) are low-luminosity X-ray 
systems (one of which has been found to be an X-ray vari- 
able AGN: J033246.8-275120; see Tables 1-2). This shows 
that the IR-excess galaxy population is heterogenous, in qual- 
itative agreement with several studies of IR-excess galax- 
ies (e.g., Teplitz etal.2007; Alexander etal.2008; Murphy 
etal.2009; Fadda etal.2010; Georgakakis etal.2010; Geor- 
gantopoulos et al. 201 1). 

3.1.3. X-ray spectral analyses 

To gain more insight into the intrinsic AGN properties (e.g., 
Nh, r, reflection components, Fe K emission) of the X-ray 
detected BzK galaxies we fitted the X-ray data using physi- 
cally motivated AGN models. Our main focus here is to con- 
strain the X-ray spectral properties of the 1 1 heavily obscured 
AGNs to identify any potential Compton-thick AGN signa- 
tures (dominant reflection component; strong Fe K emission). 
However, we also explored the X-ray spectral properties of 
all of the X-ray detected BzK galaxies to provide constraints 
on key spectral parameters (e.g., the intrinsic X-ray spectral 
slope; T) and search for further heavily obscured AGNs not 
identified using the simple X-ray spectral slope criteria. We 
extracted the X-ray spectra of each source following §2.3 and 
initially fitted each X-ray spectrum over the observed-frame 
0.5-8 keV energy band with an absorbed power-law model 
(the model components are wabs*zwabs*pow in xspec) 
using the C statistic (Cash 1979). 

We first focus on the results obtained for the nine lumi- 
nous AGNs (Lx > 10 43 erg s _1 ) since the good photon statis- 
tics (« 220^-600 net counts; mean of ss 1600 net counts) 
provide accurate constraints on the intrinsic X-ray spectral 
slope and the presence of any absorption. The individual 
best-fitting parameters for the luminous AGNs are consistent 
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FIG. 4. — Best-fit parameters (T versus Nu) from jointly fitting the X-ray 
spectra of the heavily obscured AGNs and low-luminosity X-ray systems with 
an absorbed power-law model. The contours refer to the 68%, 95% and 99% 
confidence limits and the crosses indicate the best-fitting parameters. The 
best-fitting T for the heavily obscured AGNs is inconsistent with the intrinsic 
X-ray spectral slope found for AGN activity (r Ri 1.3-2.5; see Fig. 2), sug- 
gesting that the observed X-ray emission for the overall sample is not well 
characterized by absorbed power-law emission. 



with those expected for relatively unobscured AGNs (r« 1.8 
and A^h < 10 cm -2 ). To provide tighter overall constraints 
we also performed joint spectral fitting for all of the sources, 
which determines the best-fitting F and A^h for the whole 
sample. Jointly fitting F and A^h but leaving the normalisa- 
tion of each source to vary, we obtained F = 1.71 +004 with 
low intrinsic absorption (N H = (O.65+012) x 10 22 cm -2 ); the 
best-fitting parameters obtained from % 2 fitting (r = 1 .781q q 4 

andA^H = (0.79+o Jq) x 10 22 cm -2 ) are statistically consistent 
with those obtained using the C statistic. These constraints 
can help interpret the X-ray spectral properties of the heav- 
ily obscured AGNs. The results obtained for the 27 low- 
luminosity X-ray systems are similar to those obtained for 
the luminous AGNs but with considerably larger uncertain- 
ties (« 3-180 net counts; mean of « 30 net counts); the best- 
fitting parameters from jointly fitting the X-ray spectra are 
r = 2.22^;?J with N H = (1.46^)) x 10 22 cm -2 (see Fig. 4). 

We now focus on the results obtained for the 1 1 heavily ob- 
scured AGNs. The only heavily obscured AGN that is bright 
enough for reasonable-quality individual X-ray spectral con- 
straints is J033222. 5-274603 (ss 560 net counts), which has 
r = 0.93+U; 2 ^ and N H = (0.82^) x 10 23 cm -2 ; the best- 
fitting parameters obtained from % 2 fitting are F = 1.03+q 

and A^h = (1.80j^f) x 10 23 cm -2 , with a reduced % 2 of 
1.13 for 24 degrees of freedom. The absorbing column den- 
sity of this object suggests that it is Compton thin; how- 
ever, the best-fitting X-ray spectral slope is flat. Examination 
of the residuals shows that the model significantly deviates 
from the data at observed-frame < 1.5 keV. Fitting the X-ray 
data between observed-frame 1 .5-8 keV gives F = 1 .40+qS} 
and N H = {2.05±°{%[) x 10 23 cm -2 ; the best-fitting param- 
eters obtained from % 2 fitting are F — 1 .59+o3g and Nh = 
(2.68+};^) x 10 23 cm -2 , with a reduced % 2 of 1.03 for 20 de- 
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FIG. 5. — Flux density versus rest-frame energy showing the composite rest-frame 2-20 keV spectra for the heavily obscured AGNs (filled circles; all objects 
except the X-ray bright AGN J033222. 5—274603) and low-luminosity X-ray systems (open triangles) as compared to an unabsorbed power-law model (dotted 
line; T = 2.0), a pure reflection model (solid curve; T = 1.7), and the best-fitting model to the reflection-dominated z~ AGN Swift J0601. 9-8636 (dot-dashed 
curve; Ueda etal. 2007); see §3.1.4. The inset panel shows the stacked X-ray spectra at rest-frame 2-10 keV of the heavily obscured AGNs with spectroscopic 
redshifts (filled squares: all objects; open circles: all objects except J03 3 23 5. 7-2749 16, which has been individually identified with Fe Ka emission; Feruglio 
et al. 201 1); the dashed line shows the expected rest-frame energy of Fe Ka. The properties of the heavily obscured AGNs are consistent with those expected for 
reflection-dominated systems and ss 10-50% are likely to be Compton-thick AGNs (see §3.1.4). 



grees of freedom. These best-fitting parameters now provide a 
better characterisation of the data and are consistent with that 
expected for an obscured Compton-thin AGN. The X-ray- 
8 /jm luminosity ratio of w 0.02 is also consistent with that 
expected for a Compton-thin AGN and is an order of mag- 
nitude larger than the luminosity ratio for the other heavily 
obscured AGNs (see Fig. 3). 

The best-fitting parameters for the other heavily obscured 
AGNs (« 20-120 net counts; mean of m 60 net counts) 
are determined from jointly fitting the X-ray spectra. Us- 
ing the absorbed power-law model, the best-fitting parame- 
ters are T = 0.35+^ with /V H = (5.79+^ 4 ) x 10 22 cm" 2 ; 
see Fig. 4. Such a flat intrinsic X-ray spectral slope is 
inconsistent with that found for typical AGNs and is also 
inconsistent with that expected from the inverse Compton 
scattering of accretion-disk photons (i.e., the X-ray emitting 
"corona"; Haardt & Maraschi 1993; Mushotzky, Done, & 
Pounds 1993; Reynolds & Nowak 2003). However, these 
properties are consistent with a reflection-dominated spec- 
trum, such as that typically identified in Compton-thick AGNs 
and some heavily obscured Compton-thin AGNs (e.g., George 
& Fabian 1991; Matt etal. 1996, 2000; Ueda etal. 2007; 
Eguchi et al. 2009; Comastri et al. 2010). Indeed, from jointly 



fitting the X-ray spectra of the heavily obscured AGNs with a 
reflection-dominated spectrum (the PEXRAV model in XSPEC; 
Magdziarz & Zdziarski 1995), leaving the intrinsic X-ray 
spectral slope as the jointly fitted parameter but allowing the 
normalisation of each source to vary, we obtain T= 1 .69^ ^'-, 
the reflection parameter in PEXRAV is fixed to R = — 1 (to 
produce only the reflection component), the cut-off energy is 
fixed to E_cut = 128 keV (e.g., Malizia etal. 2003), and all 
of the other parameters (inclination angle, elemental abun- 
dances) are fixed at their default values. The best-fitting X-ray 
spectral slope is now in good agreement with that found for 
the luminous AGNs, providing evidence that these systems 
are the heavily obscured reflection-dominated counterparts of 
the luminous AGNs. 

3.1.4. Reflection-dominated heavily obscured AGNs 

To explore further whether the reflection-dominated model 
provides a good description of the X-ray spectra of the 
heavily obscured AGNs, we also produced a composite 
rest-frame 2-20 keV spectrum following §3.4 of Alexander 
et al. (2005a). Briefly, the unbinned spectrum of each object is 
fitted using a simple power-law model and an unfolded spec- 
trum is produced, taking into account the Chandra effective 
area and exposure time. Each spectrum is then converted to 
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rest-frame energies and all of the spectra are combined and 
binned to increase the signal-to-noise ratio. The composite 
X-ray spectrum of the heavily obscured AGNs is shown in 
Fig. 5 and is compared to the PEXRAV model with T = 1.7. 
The similarity between the composite X-ray spectrum and the 
pure reflection model is striking, directly showing that the typ- 
ical X-ray spectrum of the heavily obscured AGNs is reflec- 
tion dominated. The composite X-ray spectrum is also consis- 
tent with that of reflection-dominated AGNs, such as those re- 
cently identified at z ~ using > 10 keV observatories {Swift; 
Suzaku; e.g., Ueda et al. 2007; Eguchi et al. 2009; Comastri 
etal. 2010); see Fig. 5. The similarity between the compos- 
ite spectrum of the low-luminosity X-ray sources and a F = 2 
power-law spectrum suggests that many of these sources are 
either intrinsically weak AGNs (see §3.1.1 for X-ray variabil- 
ity constraints) or dominated by HMXBs; see Fig. 2 for the 
typical range of X-ray spectral slopes for AGNs and HMXBs. 

If the X-ray emission of the heavily obscured AGNs is dom- 
inated by reflection then we would also expect to identify 
strong Fe K emission (e.g., Reynolds & Nowak 2003). There 
is no clear evidence for Fe K emission in the stacked spectrum 
but since half of heavily obscured AGNs have photometric 
redshifts, the Fe K emission may be smeared out. To test this 
hypothesis we only stacked the X-ray spectra of the heavily 
obscured AGNs with spectroscopic redshifts; see inset panel 
in Fig. 5. Encouragingly, we now identify a strong emission 
feature (« 1 keV rest-frame equivalent width) at « 6.4 keV, 
which is likely to be due to Fe K emission and suggests the 
presence of Compton-thick AGNs, which typically have Fe K 
emission with an equivalent width of £ 1 keV (e.g., George 
& Fabian 1991; Matt et al. 1996, 2000; Delia Ceca et al. 
2008). However, this feature is weaker (« 0.5 keV rest-frame 
equivalent width) when we remove J033235. 7-2749 16, which 
has been individually identified with strong Fe K emission 
(Feruglio et al. 201 1). Since Compton-thin AGNs have lower 
equivalent width Fe K emission than Compton-thick AGNs 
( < 0.5 keV; Mushotzky, Done, & Pounds 1993; Risaliti 
et al. 2002; Dadina et al. 2008), the weaker Fe K emis- 
sion with J033235.7-274916 removed suggests that the heav- 
ily obscured AGNs comprise a combination of reflection- 
dominated Compton-thick and Compton-thin AGNs (e.g., 
Matt et al. 2000; Ueda et al. 2007; Eguchi et al. 2009; Comas- 
tri et al. 2010). Under the assumption that the Compton-thick 
AGNs have Fe K equivalent widths of £ 1 keV, this sug- 
gests that ^5 50% of the heavily obscured AGNs are absorbed 
by Compton-thick material; conversely, a lower limit to the 
Compton-thick AGN fraction is « 10% due to the identifi- 
cation of J033235.7-274916 (Feruglio et al. 201 1). Qualita- 
tively similar results have been obtained by Georgakakis et al. 
(2010) for IR-excess galaxies at z ~ 1. 

Accurate measurements of the intrinsic X-ray luminosity 
are difficult in reflection-dominated AGNs due to large uncer- 
tainties on the reflecting geometry. We therefore employ sev- 
eral different approaches that should bracket the likely range 
in intrinsic X-ray luminosities. Clearly, a lower limit is ob- 
tained from the observed 2-8 keV luminosities, which gives 
an average luminosity of log(Lx/erg s _1 )« 43.1 at the av- 
erage rest-frame energy of w 6-24 keV. Conversely, an up- 
per limit is obtained by assuming that the rest-frame 8 fim 
emission is dominated by AGN activity (i.e., on the basis 
of the AGN-dominated line in Fig. 3); following this ap- 
proach we obtain an average intrinsic X-ray luminosity of 
log(L 2 -iokev/erg s _1 )« 44.0 for the median rest-frame 8 /xm 
luminosity of log(Ls ^/Lq) ss 11.1. Lastly, we can esti- 



mate the intrinsic X-ray luminosity of the heavily obscured 
AGNs under the reasonable assumption that they are the ab- 
sorbed counterparts of the luminous AGNs (see §3.1.3 for 
some evidence). Using the average X-ray-8 luminos- 
ity ratio of the luminous AGNs (« 0.07; a factor « 3 less 
than the intrinsic ratio of « 0.21), we estimate an aver- 
age intrinsic luminosity for the heavily obscured AGNs of 
log(L2-iokev/ er g s _1 )« 43.5. On the basis of this approach 
we predict a rest-frame 2-10 keV luminosity within a factor 
« 2 of that measured from the optical and X-ray spectroscopy 
for the Compton-thick AGN J033235.7-274916 (Feruglio 
etal. 2011; see Fig. 3), one of the heavily obscured AGNs in 
our sample. Since the average X-ray-8 pm luminosity ratio 
of the luminous AGNs is lower than that expected for AGN- 
dominated systems, a natural prediction of this final approach 
is that the rest-frame 8 pm emission from many of the heav- 
ily obscured AGNs is dominated by star-formation activity, 
in agreement with that found from Spitzer-IRS spectroscopy 
(e.g., Murphy et al. 2009; Donley et al. 2010; Fadda et al. 
2010). 

3.2. X-ray undetected BzK galaxies 

The majority of the BzK galaxy population remains unde- 
tected in the 4 Ms Chandra observation (see §3.1) but may 
host X-ray weak AGN activity below the detection limit. We 
can place constraints on the presence of heavily obscured 
AGN activity in these systems using X-ray stacking analyses. 
We stacked the X-ray data of the X-ray undetected BzK galax- 
ies adopting the procedure in §2.5. The X-ray stacking results 
are presented in Table 3. As previously found from the X- 
ray stacking analyses of BzK galaxies in Daddi et al. (2007a), 
the result for the IR-excess galaxies differs from that of the 
IR-normal galaxies: Y = 1.4+Q3 for the IR-excess galaxies 

and r = 2.0^|j 4 for the IR-normal galaxies; see Table 3. The 
stacked X-ray spectral slope for the IR-normal galaxies is con- 
sistent with that found by Daddi et al. (2007a) using the 1 Ms 
Chandra data (r « 1.8) but the X-ray spectral slope for the 
IR-excess galaxies is significantly steeper (r w 0.9 was ob- 
tained by Daddi et al. 2007a). 

The lack of a flat X-ray spectral slope for the IR-excess 
galaxies appears to suggest that we have now individually de- 
tected many of the heavily obscured AGNs that were origi- 
nally contributing to the stacked data in Daddi etal. (2007a). 
However, we must be careful when interpreting X-ray stack- 
ing analyses of X-ray undetected source populations since the 
effects of Eddington bias and source variability (see §3.1.1) 
can dominate over the X-ray signal produced by the major- 
ity of the source population. For example, in going from 
the 1 Ms Chandra data to the 4 Ms Chandra data, the same 
number of heavily obscured AGNs were identified in the IR- 
excess (three of the 13 X-ray detected sources) and IR-normal 
(three of the 14 X-ray detected sources) galaxy populations 
despite there being little evidence for heavily obscured AGNs 
in the IR-normal galaxies from the 1 Ms stacked data of Daddi 
etal. (2007a). It is therefore likely that further heavily ob- 
scured AGNs will be identified with deeper Chandra data but 
we cannot provide direct X-ray constraints from the current 
dataset. 

3.3. Re-evaluation of the space density of distant 
Compton-thick AGNs 

Distant Compton-thick AGNs are of great scientific in- 
terest since they may produce a large fraction of the unre- 
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FIG. 6. — Space density of heavily obscured AGNs with L\ i 10 43 ergs -1 . 
The plotted data only corresponds to the results derived from the X-ray 
detected heavily obscured AGNs identified here (filled squares), candidate 
Compton-thick AGNs identified from X-ray spectral analyses in the CDF- 
S (filled triangles; Tozzi et al. 2006), and X-ray stacking analysis results 
of X-ray undetected candidate Compton-thick AGNs with Lx i 10 43 erg s 
(open squares; Daddi et al. 2007a; Fiore et al. 2008; Treister et al. 2009). The 
solid bar indicates the space-density estimates for a range of Compton-thick 
AGN percentages for the heavily obscured AGNs identified here. These re- 
sults are compared to the space-density predictions for Compton-thick AGNs 
with Lx ^ 10 43 erg s _I based on the models of Gilli et al. (2007; dashed 
curve) and Triester, Urry, & Virani (2009; dotted curve). 



solved > 8 keV background (e.g., Worsley et al. 2005; Gilli 
et al. 2007; Treister, Urry, & Virani 2009). In §3.1 we 
used a variety of analyses (strong reflected-dominated spec- 
trum, identification of Fe K, small X-ray-8 luminosity 
ratios) to infer that « 10-50% of the X-ray detected heav- 
ily obscured AGNs are Compton thick; the lower limit on the 
Compton-thick AGN fraction corresponds to the clear iden- 
tification of the Compton-thick AGN in J033235.7-274916 
(Feruglio et al. 201 1) while the upper limit corresponds to the 
constraints derived from the stacked X-ray spectrum when 
J033235. 7-274916 is removed. We can use these constraints 
to better estimate the space density of distant Compton-thick 
AGNs. In this calculation we have assumed a broad red- 
shift range of z — 1.4-2.6, which gives a comoving volume 
of « 0.7 Gpc 3 for the 5.5 arcmin radius region explored here. 
We have taken into account the 10% incompleteness in the 
BzK galaxy selection due to blended Spitzer-IRAC sources 
and the 30% incompleteness due to unreliable UV slopes (see 
§5.1 in Daddi etal. 2007a). 

Following the procedure outlined above, we calculate a 
Compton-thick AGN space density of Oc-thick ~ / x 4 x 
10~ 5 Mpc~ 3 at z ~ 1.4-2.6, where f corresponds to the 
Compton-thick AGN fraction in our heavily obscured AGN 
sample; see Fig. 6. As argued in §3.1.4, the average in- 
trinsic X-ray luminosity of these heavily obscured AGNs is 
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constraints derived from X-ray stacking analyses of X-ray 
undetected IR-bright galaxy populations (Daddi etal. 2007a; 
Fiore etal. 2008): our space density estimates are >, 10 
times lower than those of Daddi etal. (2007a), who used the 
same object-selection approach as that adopted here but relied 
only on X-ray stacking analyses of X-ray undetected galax- 
ies. Our space-density constraints are also broadly consistent 
with those of Tozzi et al. (2006), who identified reflection- 
dominated AGNs using X-ray spectral analyses of sources 
detected in the 1 Ms CDF-S observations; we have plot- 
ted the Tozzi et al. (2006) space density in Fig. 5 assuming 
that 100% of the reflection-dominated AGNs are Compton 
thick and therefore this represents the maximum space den- 
sity from that study. However, our space-density estimate 
is also a lower limit on the true Compton-thick AGN space 
density since (1) there may be further Compton-thick AGNs 
with intrinsic i2-i0kev ~ 10 43 erg s _1 that lie below our X- 
ray detection limit (i.e., those with A/h S> 10 24 cm~ 2 and a 
weak reflection component; e.g., Matt et al. 2000), and (2) our 
sample does not include distant Compton-thick AGNs not se- 
lected using the BzK technique. An example of the latter is the 
z = 1.53 X-ray bright AGN CXO J033218. 3-275055, which 
has strong Fe K emission identified in the 3 Ms XMM-Newton 
observations of the CDF-S (Comastri etal. 2011) but is not 
selected as a BzK galaxy. 5 A more complete AGN selection 
can be derived using mid-to-far-infrared selection, which is 
the focus of a future paper (A. Del Moro et al. in prep). 

Many studies have predicted the space density of distant 
Compton-thick AGNs from the X-ray luminosity functions 
of relatively unobscured AGNs and X-ray background con- 
straints (see Ballantyne et al. 201 1 for a comparison of many 
of the current studies). In Fig. 6 we compare our space- 
density constraints with the predictions from Gilli et al. (2007) 
and Treister, Urry, & Virani (2009) for L x > 10 43 erg s -1 , 
which broadly represent the most optimistic and pessimistic 
estimates, respectively. The model predictions are already in 
broad agreement with our range of space-density measure- 
ments, despite our conservative source-selection approach. 
However, given the significant uncertainties in the fraction 
of Compton-thick AGNs and sample incompleteness, strong 
conclusions cannot be derived from the current data. 

We can also compare our derived Compton-thick AGN 
space density to constraints for other distant AGN popu- 
lations. From a variety of studies, the measured space 
density of Lx ^ 10 43 erg s _1 AGNs at z ~ 2 ranges 
from « (1-2) xl0~ 5 Mpc~ 3 for unobscured AGNs (Nu < 
10 22 cm" 2 ) to w (1-7) xl0~ 5 Mpc~ 3 for all AGNs (e.g., 
Barger et al. 2005; Hasinger et al. 2005; La Franca et al. 2005 ; 
Silverman etal. 2008; Aird etal. 2010). On the basis of 
these studies, our constraints imply that the space density of 
Compton-thick AGNs at z ~ 2 is comparable to that of unob- 
scured AGNs at z ~ 2 and also suggests that Compton-thick 
AGNs comprise a non-negligible fraction of the AGN popula- 
tion at z 2. However, our constraints do not yet support the 
hypothesis that Compton-thick AGNs outnumber Compton- 
thin AGNs at high redshift. The direct identification of in- 
dividual Compton-thick AGN signatures from X-ray, optical, 
and mid-IR spectroscopy (e.g., Alexander et al. 2008; Comas- 
tri etal. 2011; Feruglio etal. 2011; Goulding etal. 2011) in 
a statistically significant number of objects ( £ 10-20 ob- 

5 We also note as an aside that only five (s=s 45%) of the 11 X-ray de- 
tected heavily obscured AGNs would be selected using the Fiore et al. (2008) 



our most optimistic space-density estimates lie below the selection criteria of R-K> 4.5 and /Mum//* > 1000. 
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jects) are required to provide more reliable constraints. New 
facilities such as NuSTAR (high-energy 6-78 keV imaging; 
Harrison et al. 20 1 0) and JWST (optical-mid-IR spectroscopy; 
Gardiner etal. 2006) may also provide improved constraints, 
along with deeper ChandralXMM-Newton observations and 
future proposed X-ray observatories such as WFXT (Murray 
etal. 2010). 

4. SUMMARY 

We have used the 4 Ms Chandra Deep Field-South obser- 
vation to constrain the ubiquity of heavily obscured AGNs in 
the z ~ 2 BzK galaxy population. Our main results are: 

• Forty seven of the 222 BzK galaxies are X-ray de- 
tected in the central region of the 4 Ms CDF-S field: 
1 1 are heavily obscured AGNs (r 1), 9 are luminous 
AGNs (Lx ~ 10 43 erg s _1 ), and 27 are low-luminosity 
X-ray systems (relatively unobscured AGNs and star- 
burst galaxies). Thirteen (w 48%) of the 27 X-ray de- 
tected BzK galaxies with reasonable-quality X-ray data 
(> 20 counts in the 0.5-8 keV band) are found to be 
variable in the X-ray band, including 6 luminous AGNs, 
4 heavily obscured AGNs, and 3 low-luminosity X-ray 
sources. See §2.2, §2.4, and §3.1.1. 

• The overall X-ray spectra of the heavily obscured 
AGNs are better characterized by a pure reflection 
model than an absorbed power-law model, suggesting 
extreme Compton-thick absorption (Na ~ 10 24 cm -2 ) 
in many systems. The identification of an emission- 
line feature at rest-frame w 6.4 keV in the composite 
2-20 keV spectrum and the small X-ray-8 /jm lumi- 
nosity ratios for the majority of these systems provide 
further support for this interpretation. See §3.1 .2-3.1 .4. 

• Many of the heavily obscured AGNs are IR-excess 
galaxies. However, only w 25% of the X-ray detected 
IR-excess galaxies are heavily obscured AGNs, which 
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is otherwise composed of relatively unobscured AGNs 
and starburst galaxies. See §3.1.2. 

• X-ray stacking analyses of the X-ray undetected BzK 
galaxies do not clearly reveal the presence of further 
X-ray undetected AGNs below the Chandra detection 
limit. This does not rule out the possibility that many 
other heavily obscured AGNs will be detected with 
deeper X-ray observations but it does suggest that they 
are not the dominant X-ray undetected population. See 
§2.5 and §3.2. 

• We estimate a Compton-thick AGN space density of 
*c-thick ~/x4x 10~ 5 Mpc~ 3 at z ~ 1.4-2.6, where 
/ lies between « 0.1-0.5. Although highly uncertain, 
these constraints are already consistent with the range 
of predictions from X-ray background models and im- 
ply that the space density of Compton-thick AGNs at 
Z ~ 2 is comparable to that of unobscured AGNs at 
z«2. See §3.1.4 and §3.3. 
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Table 1 . Overall Properties of the X-ray detected BzK Galaxies ^ 

to 



A'-band X-K Xyega ^uv ^2-lOkeV vL S)im SFR 1 Ms 

ay2000 a 8j2000 a (arcsec)* f (mag)" XID d log(L. 3 ) a log(erg s" 1 / log(L e ) fl Excess" F* Source?' Notes^ 



03 32 10.95 


—27 48 56.1 


0.48 


2.81 


20.2 


167 


12.04 


43.96 


11.72 


17.3 


03 32 12.55 


-27 49 38.2 


0.29 


2.45 


21.3 


185 


11.54 


42.14 


10.64 


1.3 


03 32 14.12 


-27 49 10.2 


0.59 


2.18 


20.6 


202 


11.33 


42.06 


11.20 


15.0 


03 32 14.42 


-27 51 10.7 


0.39 


1.544 


19.2 


205 


12.11 


41.98 


10.72 


0.5 


03 32 14.79 


-27 44 02.5 


0.67 


1.56 


20.7 


208 


10.80 


42.33 


10.37 


2.9 


03 32 16.94 


-27 50 04.0 


0.35 


1.613 


20.6 


236 


11.40 


41.59 


10.50 


1.1 


03 32 17.81 


-27 52 10.3 


0.57 


1.76 


21.7 


247 


10.98 


<41.95 


10.47 


2.7 


03 32 18.24 


-27 52 41.2 


0.45 


2.801 


21.4 


254 


11.35 


43.41 


10.74 


2.9 


03 32 21.30 


-27 51 01.5 


0.52 


1.84 


20.2 


293 


11.70 


41.68 


11.10 


4.5 


03 32 21 99 


-27 51 11.9 


0.28 


3.64 


21.2 


298 


11.39 


42.76 


11.85 


124.6 


03 32 22.54 


-27 46 03.8 


0.44 


1.730 


20.4 


308 


11.26 


42.90 


11.00 


8.7 


03 32 22.55 


-27 48 14.9 


0.48 


2.54 


20.4 


310 


11.63 


41.90 


11.20 


7.4 


03 32 24.84 


-27 50 50.1 


0.38 


2.28 


20.7 


337 


11.47 


41.91 


11.14 


8.7 


03 32 25.98 


-27 47 51.3 


0.69 


1.90i 


20.5 


360 


11.47 


41.91 


11.12 


8.0 


03 32 28.79 


-27 47 55.5 


0.95 


1.383 


19.3 


394 


12.07 


41.59 


10.75 


0.6 


03 32 29.09 


-27 46 29.0 


0.34 


2.227 


19.8 


399 


11.89 


41.91 


11.27 


5.2 


03 32 29.48 


-27 43 22.0 


0.56 


1.609 


19.8 


405 


11.66 


<41.90 


10.66 


1.1 


03 32 29.99 


-27 45 29.9 


0.38 


1.218 


18.3 


417 


11.69 


43.74 


11.03 


3.6 


03 32 31.47 


-27 46 23.2 


0.43 


2.223 


19.0 


435 


12.41 


42.43 


11.71 


7.2 


03 32 31.52 


-27 48 53.8 


0.43 


1.879 


20.3 


437 


11.66 


42.41 


11.25 


8.3 



1 fiS+ 008 
1 - bS -0.0£) 

>0.53 
>0.57 

1 3+°- 23 
u - ,: >-0.20 

l -^ J -0.29 

> 0.60 

1 -7Q+0.26 

1 -' y -0.20 

36+ 64 
u - ,D -0.44 
om +0.21 

17+ 007 

-0 S9+ - 26 
uoy„ 22 

>0.43 

>0.94 

> 0.64 

> 0.42 
<0.09 

1 44+0.04 
lM -0.04 

1 00+ 032 
l.UU_Q 27 

1 31+0- 34 
1 --' 1 _0.27 



Y 
N 
N 
Y 
N 

N 
N 
Y 
N 
Y 

Y 

N 
N 
N 
N 

N 
N 
Y 
Y 
Y 



03 32 31.55 


-27 50 28.6 


0.76 


1.613 


18.9 


436 


12.14 


<41.64 


11.01 


1.2 


<0.70 


N 


03 32 32.93 


-27 50 40.5 


0.82 


2.50 


20.5 


451 


11.96 


41.92 


11.08 


2.3 


> 0.23 


N 


03 32 34.04 


-27 50 28.7 


0.44 


1.384 


19.5 


467 


11.76 


41.42 


10.71 


1.0 


>0.46 


N 


03 32 34.46 


-27 50 04.9 


0.21 


2.14 


20.5 


474 


11.46 


41.99 


11.41 


22.5 


> 0.77 


N 


03 32 34.98 


-27 49 31.9 


0.97 


2.55 


21.9 


481 


11.44 


42.03 


11.30 


16.1 


> 0.51 


Y 



Heavily obscured AGN; optical AGN 



Optical AGN 

Heavily obscured AGN 

Heavily obscured AGN 

Heavily obscured AGN 
Heavily obscured AGN 

IRS: s/burst dominated 



Heavily obscured AGN 
Optical AGN 
Heavily obscured AGN 
IRS: s/burst dominated 

Heavily obscured AGN 



03 32 35.72 


-27 49 16.1 


0.31 


2.578 


20.0 


490 


11.85 


42.51 


12.25 


171.0 


03 32 35.97 


-27 48 50.4 


0.31 


1.309 


19.1 


493 


11.93 


42.24 


10.55 


0.4 


03 32 36.17 


-27 51 26.5 


0.37 


1.613 


20.0 


499 


11.69 


42.67 


10.80 


1.6 


03 32 36.18 


-27 46 27.6 


0.43 


2.48 


20.9 


501 


11.96 


41.99 


11.03 


1.9 


03 32 37.36 


-27 46 45.5 


0.62 


1.843 


20.0 


512 


11.71 


41.74 


11.00 


3.0 


03 32 37.74 


-27 50 00.6 


0.81 


1.619 


19.2 


517 


12.23 


41.71 


11.16 


1.6 


03 32 37.77 


-27 52 12.3 


0.27 


1.603 


18.8 


518 


11.98 


44.16 


11.64 


15.1 


03 32 37.96 


-27 53 07.9 


0.26 


1.97 


21.2 


520 


11.44 


42.60 


10.55 


1.2 


03 32 38.55 


-27 46 34.2 


0.77 


2.55i 


21.3 


525 


11.68 


42.27 


11.77 


46.7 


03 32 39.08 


-27 46 02.1 


0.17 


1.216 


18.3 


537 


11.85 


43.22 


10.82 


1.2 


03 32 39.74 


-27 46 11.5 


0.08 


1.552 


19.4 


549 


11.63 


43.24 


11.08 


4.9 


03 32 40.06 


-27 47 55.4 


0.45 


1.998 


19.5 


552 


12.08 


41.99 


11.58 


9.8 


03 32 40.76 


-27 49 26.2 


0.40 


2.130 


19.7 


555 


11.85 


42.05 


11.34 


7.2 


03 32 41.80 


-27 51 35.3 


0.20 


1.63 


20.2 


562 


11.54 


41.90 


10.94 


3.7 


03 32 43.25 


-27 49 14.3 


0.20 


1.920 


19.5 


577 


11.31 


43.96 


11.10 


10.9 


03 32 43.46 


-27 49 01.8 


0.46 


1.78i 


20.2 


579 


11.76 


42.00 


11.27 


7.0 


03 32 43.61 


-27 46 59.0 


0.20 


1.570 


20.4 


580 


11.46 


41.86 


10.66 


1.7 



n 44+0.23 
-0.20 
1 98+ 043 
, ,,7+0.18 

>0.41 
>0.42 



> 0.64 

, -7 1+0.03 
1 -' 1 -0.03 

> 1.41 

> 0.86 

i 1^+0.05 
1 - 1J -0.05 



S+0.07 
-0.07 

>0.63 
>0.65 

77+0- 57 
"•"-0.38 

1 cc+0.04 
1 - JJ -0.04 



> 0.88 

> 0.83 



Y 
Y 
Y 

N 
N 

N 
Y 
Y 
N 
Y 

Y 
N 
N 
N 
Y 

Y 
N 



Heavily obscured AGN; optical AGN 



Optical AGN; IRS: AGN dominated 

IRS: s/burst dominated 
Optical AGN 

IRS: s/burst dominated 
IRS: s/burst dominated 
IRS: s/burst dominated 
Heavily obscured AGN 
Optical AGN 

IRS: s/burst dominated 



Table 1. Overall Properties of the X-ray detected BzK Galaxies — Con- 
tinued 



K-bwd 

0/2000" 8/2000" 


X-K 
(arcsec)* 


z c 


-Kvega 

(mag)" 


XID"* 


log(Z, Q )° 


-f-2-lOkeV 
log(erg s _1 ) d 


log(L ) a 


SFR 

Excess" 




1 Ms 
Source? 15 


Notes^ 


03 32 44.02 
03 32 44.37 
03 32 44.60 


-27 46 34.9 
-27 49 11.3 
-27 48 36.0 


1.00 
0.53 
0.18 


2.688 
2.13 
2.593 


20.9 
20.4 
21.4 


583 
589 
593 


11.75 
11.25 
11.55 


43.63 
<42.03 
43.38 


11.69 
11.44 
11.11 


31.1 
40.5 
6.6 


2 00+ 015 
z - uu -0.16 

>0.60 

1 in+0.12 

L - ly -o.u 


Y 
N 
Y 


IRS: s/burst dominated 


03 32 46.84 
03 32 47.72 


-27 51 20.9 
-27 50 38.0 


0.25 
0.40 


2.292 
1.63 


20.4 
19.2 


617 
625 


11.74 
12.06 


42.42 
41.94 


11.06 
11.21 


3.6 
2.8 


> 0.73 

33+ - 20 


N 
N 


Heavily obscured AGN 



a BzK galaxy properties taken from Daddi etal. (2007b). Co-ordinates correspond to the /f-band position of the BzK galaxy. The UV luminosity corresponds to rest-frame 1500 A and has been corrected for 
extinction (see §3.6 of Daddi et al. 2007b). The rest-frame 8 fjm luminosity is calculated using the 24 fjm flux density, with small K-corrections applied (see §3. 1 of Daddi et al. 2007b). The SFR excess corresponds to 
the ratio of star-formation rates (mid-IR+UV versus extinction-corrected UV; see §2.2 of Daddi etal. 2007a). 

b Offset between the position of the X-ray source and the /f-band BzK counterpart in arcseconds. 

c BzK galaxy redshifts. Optical spectroscopic redshifts are denoted by having three decimal places and come from Szokoly et al. (2004), Mignoli et al. (2005), Cimatti et al. (2008), Vanzella et al. (2008), Popesso 
et al. (2009), and Balestra et al. (2010). All other redshifts are either photometric (from Grazian etal. 2006; Luo et al. 2010; Rafferty et al. 2011) or from Spitzer-IRS spectroscopy (highlighted with "i" and from 
Teplitz et al. 2007 and Fadda et al. 2010). 

d X-ray source properties. XID corresponds to the X-ray identification number in Xue et al. (2011). The X-ray spectral slope (r) determined from the band ratio (2-8 keV to 0.5-2 keV count-rate ratio) plus 
1 c uncertainty are taken from Xue etal. (2011); T is re-calculated from the band ratio for the low-count sources listed in Xue etal. (2011) with T = 1.4. The rest-frame 2-10 keV luminosity is calculated from the 
observed-frame 0.5-2 keV flux in Xue et al. (201 1) and converted to rest-frame 2-10 keV assuming T = 1.8. 

e Indicates if the source was detected in the 1 Ms CDF-S catalogs of Alexander et al. (2003); otherwise the source was detected in the 4 Ms CDF-S catalogs of Xue et al. (201 1). 

f Notes and classifications. Objects classified as "Heavily obscured AGN" have r ,$ 1 and objects classified as optical AGN have optical spectroscopic signatures from Szokoly et al. (2004) indicating AGN activity. 
We have also indicated which sources have Spitzer-IRS spectroscopy, available either from Teplitz et al. (2007), Donley et al. (2010), or Fadda et al. (2010). 
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Table 2. X-ray Properties of the X-ray detected BzK Galaxies 



X 


ray 




X-ray variability constraints c 


Variable 


f^./2000 


OJ2000 


Net rnnnts' 3 


X 2 




a 2 

^rms 


°HMXB 


AGN? d 


03 32 10.98 


-27 48 56.5 


844 


4.09 


0.004 


0.014 


0.0054 


Y 


03 32 12.57 


-27 49 38.2 


9 












m Tl 1 A l C 

03 32 14.15 


—11 49 10.6 


1 8 












03 32 14.43 


-27 51 11.0 


82 


4.30 


0.005 


0.141 


0.0011 


Y 


03 32 14.80 


-27 44 03.2 


89 


0.371 


0.723 


-0.034 


N/A 


N 


03 32 16.96 


-27 50 04.2 


5 












03 32 17.84 


-27 52 10.8 


27 


0.952 


0.289 


-0.001 


0.0030 


N 


m io 10 O/l 

03 32 18.24 


—27 52 41.6 


215 


1.46 


0.183 


0.009 


0.0019 


N 


03 32 21.31 


-27 51 02.0 


22 


0.688 


0.416 


-0.092 


0.0110 


N 


03 32 22.00 


-27 51 12.1 


80 


0.106 


0.952 


-0.058 


0.0110 


N 


03 32 22.56 


-27 46 04.2 


564 


6.47 


0.000 


0.036 


0.0005 


Y 


03 32 22.58 


-27 48 15.2 


58 


0.228 


0.842 


-0.073 


0.0110 


N 


ni n ia oc 
Uj 51 24. 5J 


0-7 cn cn a 
— 1 1 00 3U.4 


7 












03 32 26.01 


-27 47 51.8 


8 












03 32 28.85 


-27 47 56.0 


17 












03 32 29.11 


-27 46 29.3 


12 












03 32 29.50 


-27 43 22.5 


54 


0.682 


0.451 


-0.038 


N/A 


N 


i\i n an c\r\ 

03 32 29.99 


—27 45 30.3 


3916 


130 


0.000 


0.163 


0.0003 


Y 


03 32 31.48 


-27 46 23.6 


59 


0.870 


0.352 


-0.017 


0.0000 


N 


03 32 31.55 


-27 48 54.0 


61 


0.670 


0.481 


-0.046 


0.0081 


N 


03 32 31.51 


-27 50 29.0 


16 












03 32 32.95 


-27 50 41.3 


10 












03 32 34.03 


CA on 1 

-27 50 29. 1 


12 












03 32 34.47 


-27 50 05.1 


10 












03 32 35.04 


-27 49 32.5 


4 












03 32 35.72 


-27 49 16.4 


70 


1.64 


0.144 


0.019 


0.0110 


N 


03 32 35.98 


-27 48 50.7 


84 


5.26 


0.001 


0.300 


0.0014 


Y 


fY2 n 16. 1 O 

\)5 51 jo. lo 


— 11 M zo.o 


179 


2.13 


0.067 


0.027 


0.0009 


N 


03 32 36.19 


-27 46 28.0 


3 












03 32 37.38 


-27 46 46.1 


11 












03 32 37.75 


-27 50 01.4 


14 












03 32 37.77 


-27 52 12.6 


4600 


250 


0.000 


0.021 


0.0001 


Y 


fV2 n n at^ 

uj J2 5 /.yo 


T7 ci no 1 
—11 jj 08.2 


71 


3.38 


0.014 


0.141 


0.0004 


Y 


03 32 38.52 


-27 46 34.9 


14 












03 32 39.09 


-27 46 02.1 


1331 


134 


0.000 


0.346 


0.0005 


Y 


03 32 39.74 


-27 46 11.5 


760 


0.833 


0.453 


-0.001 


0.0003 


N 


03 32 40.05 


-27 47 55.8 


26 


1.50 


0.145 


-0.014 


0.0110 


N 


/"v-> -> o j j"\ n 

03 32 40.77 


-11 49 26.6 


13 












03 32 41.82 


-27 51 35.4 


37 


4.09 


0.007 


0.508 


0.0007 


Y 


03 32 43.24 


-27 49 14.5 


2117 


39.2 


0.000 


0.077 


0.0004 


Y 


03 32 43.45 


-27 49 02.2 


19 












03 32 43.61 


-27 46 59.2 


16 












03 32 44.04 


-27 46 35.9 


380 


1.13 


0.281 


0.001 


N/A 


N 


03 32 44.41 


-27 49 11.3 


7 












03 32 44.61 


-27 48 36.2 


333 


3.40 


0.011 


0.033 


0.0003 


Y 


03 32 46.86 


-27 51 21.0 


43 


4.34 


0.006 


0.335 


0.0056 


Y 


03 32 47.73 


-27 50 38.4 


118 


10.4 


0.000 


0.210 


0.0002 


Y 



a Co-ordinates correspond to the X-ray position from Xue et al. (201 1). 

b Background-subtracted (net) counts in the 0.5-8.0 keV band used in the X-ray spectral analyses and X-ray variability analyses. 

c X-ray variability constraints for the sources with reasonable-quality X-ray data (> 20 counts in the 0.5-8 keV band). The % 2 statistic is calculated by 
comparing the variability observed between observations to that expected from Poisson statistics. P^i gives the probability that the observed variability is due to 

Poisson noise, a 2 ms gives the normalized excess variance (Nandra et al. 1997), and Ohmxb gi ves me upper limit to the HMXB contribution to the variability. N/A 
indicates that no SFR information was available and so no attempt was made to estimate o" HMX b : this does not affect any significantly variable sources. See §2.4. 
d Indicates if the source was found to show excess X-ray variability over that expected from the HMXB population; see §2.4. 
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Table 3. X-ray Stacking Analyses of the X-ray Undetected BzK Galaxies 













0.5-2 keV 


2-8 keV 






^2-lOkeV 


Sample 


N 


z a 


log(Z- ) a 


log(L ) a 


(10- 6 cts s-')* 


(Kr 6 cts s- 1 )* 


Band ratio c 




log(erg s"') d 


IR excess galaxies: K < 22 


47 


1.97 ±0.34 


11.45 ±0.46 


11.01 ±0.33 


1.47 (12.6) 


0.84 (3.9) 


0.57±0.15 


1.4±0.3 


41.40 


IR normal galaxies: K < 22 


116 


1.93 ±0.44 


11.48 ±0.26 


10.53 ±0.30 


1.19 (15.8) 


0.42 (2.8) 


0.33 ±0.12 


2.0 ±0.4 


41.32 



Median galaxy properties and MAD (see Footnote 3): redshift, exinction-corrected UV luminosity (rest-frame 1500 A; see §3.6 of Daddi etal. 2007b), 
rest-frame 8 /jm luminosity (calculated using the 24 /jm flux density, with small K-corrections applied; see §3.1 of Daddi etal. 2007b). 
b Count rates in the 0.5-2 keV and 2-8 keV bands; the numbers in parantheses correspond to the S/N ratio. 

X-ray spectral properties: band ratio (2-8 keV to 0.5-2 keV count-rate ratio) and X-ray spectral slope (r), derived from the band ratio, and 1 a uncertainties. 
Calculated following Xue et al. (2011). 

d Rest-frame 2-10 keV luminosity calculated from the 0.5-2 keV flux and converted to rest-frame 2-10 keV assuming r = 1.8; the 0.5-2 keV flux is calculated 
following Xue etal. (2011). 



